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The hyperfine interaction between the quadrupole moment of atomic nuclei and the electric field 
gradient (EFG) provides information on the electronic charge distribution close to a given atomic 
site. In ferroelectric materials, the loss of inversion symmetry of the electronic charge distribution 
is necessary for the appearance of the electric polarization. We present first-principles density 
functional theory calculations of ferroelectrics such as BaTiOa, KNbOs, FbTiOa and other oxides 
with perovskite structures, by focusing on both EFG tensors and polarization. We analyze the 
EFG tensor properties such as orientation and correlation between components and their relation 
with electric polarization. This work supports previous studies of ferroelectric materials where a 
relation between EFG tensors and polarization was observed, which may be exploited to study the 
ferroelectric order when standard techniques to measure polarization are not easily applied. 



I. INTRODUCTION 

There is great interest in ferroelectric/multiferroic ma- 
terials nowadays due to their potential application in a 
plethora of subjects, ranging from high density mem- 
ories to magnetoelectric sensors »i"— The complexity of 
electronic phenomena at the nanoscale makes them a 
hot research topic and a fertile ground for new exper- 
imental techniques, which are able to probe point-like, 
atomic-scale properties. To this aim, the use of local 
probes as in hyperfine interactions techniques, such as 
Mossbauer effect spectroscopy, perturbed angular corre- 
lation (PAC) spectroscopy, nuclear magnetic resonance 
(NMR), and nuclear quadrupole resonance (NQR), give 
access to atomic scale information of the electronic charge 
density- through the measure of the electric field gradi- 
ents (EFGs), thus probing the phenomenology of mate- 
rials at the nanoscale. 

From a theoretical point of view, advances in modern 
density functional theory have made the calculation of 
the spontaneous electric polarization P a routine calcu- 
lation in an ab initio framework. ^"^ Recently, the cal- 
culation of EFGs has become also possible from first- 
principles fSr™ However, to the best of our knowledge, a 
theoretical study based on ab-initio density functional 
methods aiming to investigate both the P and EFGs in 
a ferroelectric material is still missing in the current lit- 
erature despite the fact that a linear correlation!^"— be- 
tween P and EFG values at a given site was shown a long 
time ago. This "correlation" would suggest that informa- 
tion obtained through the measurements of EFG tensors 
can provide indirect access to the polarization as well. 
If true, one could study a macroscopic property of the 
crystal, such as P by using local probes. In this work, 
we want to explore such a possibility. 

A nucleus with a nonspherical nuclear charge distri- 



bution possesses an electric quadrupole moment which 
leads to a hyperfine splitting for a nuclear spin / > 1 
if subjected to an EFG. The hyperfine techniques pre- 
viously mentioned can measure the quadrupole coupling 
constant, which is the interaction between the nuclear 
quadrupole moment and the EFG. The EFG, in turn, 
arises due to the Coulomb potential at the nucleus, and 
its measurement is sensitive to the surrounding electronic 
charge density. More precisely, it is defined as the sym- 
metric traceless second-rank tensor of second derivatives 
of the Goulomb potential with respect to the spatial coor- 
dinates, Vij = d'^V/{dxidxj), at the nuclear position. In 
the principal axis coordinate system, the tensor is diago- 
nal and its elements are usually ordered by the convention 
\Vzz\ > \Vyy\ > \Vxx\- Usually Vzz and the asymmetry 
parameter 77 = {Vxx — Vyy)/Vzz are used in the analysis 
of measurements. We recall that the EFG is site depen- 
dent, and its principal axes {x,y,z) may not be the same 
at every site, although they are usually along symmetry 
axes of the crystal. EFG studies are found in various 
types of materials, for example, intermetallics,^^ metal 
complexeSfi^ magnetic)^ or multiferroic compoundsi^ 

Previous studies have shown that in some ferro- 
electric materials P follows a temperature dependence 
which can be related to the EFG at specific atomic 
sites. For instance, NMR using ^■^Na in Rochelle salts 
[NaK(tartrate) -41120] showed that P and the EFG are 
linearly relatedJ^"— In 1978 Yeshurun suggested^^ that 
the EFG due to static displacements should be propor- 
tional to P^ in perovskite crystals using an empirical 
model for interpreting previous ^^Fe Mossbauer measure- 
ments in BaTiOa.^'^ Dynamical aspects were also consid- 
ered by relating the EFG to the electric susceptibility, 
and it was found that the EFG should have a critical be- 
havior when approaching Tc-^^ This peculiar feature was 
recently used in the identification of ferroelectricity with 
EFG measurements^ In this work, Pri^2;Caa;Mn03 was 



studied with the measurement of the EFG at ^^-'^™Cd 
probes implanted into the sample. An abrupt change 
was found in a short temperature interval. This was as- 
sociated to the onset of ferroelectricity, since the EFG 
should be dominated by a contribution proportional to 
the electric susceptibility at the transition, i.e., with its 
critical behavior.-?- In the same work,^'* it was also sug- 
gested that the temperature dependence of EFG tensors 
can give information on the onset of charge or orbital or- 
dering. In Ref.[23it was argued that the static part of V^z 
should have the following behavior with respect to P: ei- 
ther it is proportional to P^ in sites which have inversion 
symmetry in the paraelectric structure, or it is propor- 
tional to P, like in Rochelle salts. The quadratic relation 
was supported by experiments in NaN02,?" in PbHfOa 
not too close to Tc (where critical behavior is found), by 
PAC measurements, — 12& and by NMR measurements in 
BaTiOs.^"^ It is therefore clear from the current literature 
that P and the EFG tensor are closely related quantities. 

Our study aims to explore this relationship by calcu- 
lating both P and EFG for simple ferroelectric materials 
and studying a possible correlation between these quan- 
tities. Some of the previously mentioned results are ob- 
tained by impurity probes in the host materials. Here, 
we shall limit our studies to systems where the probes 
are natural constituents of the materials. 

This work is organized as follows. In Sec. |TT] we dis- 
cuss the computational details. In Sec. IIII Al we present 
the results and discuss the relationship between P and 
EFGs for simple tetragonal or orthorhombic (Sec. IIIIBp 
systems. We analyze the possible linear correlations be- 
tween EFG tensor components in Sec. IIII CI A study of 
the variation of Vzz (P) with the atomic numbers of dif- 
ferent materials is shown in section ITVl Finally, in Sec.lVl 
we draw our conclusions. 



II. TECHNICAL DETAILS 



TABLE L Lattice constants (in A) and atomic distortions 
used in the calculations of the perovskite compounds. 5a, 
<5s, <5oi, &01 correspond to the fractional distortions in the c 
direction, with respect to the ideal positions, of the inequiv- 
alent sites. Sites A are at the corners of the unit cell, B is 
the transition metal inside an octahedron formed by apical 
Ol and equatorial 02 oxygen atoms. 



Compound a 


c 5a <5s 5oi 5o2 


BaTiOa^ 4.00 


4.03 0.02 -0.03 -0.01 


KNbOa '''" 4.00 


4.06 0.02 0.04 0.04 


PbTiOs "^'^ 3.90 


4.16 0.04 0.11 0.12 


BaZrOs' 4.19 




CaTiOs' 3.83 




NaNbOs* 3.94 




PbZrOa' 4.13 




SrTiOa' 3.91 




LiNbOa'^ 4.00 




^ Reference 


2i- 




'^ Reference 


34. 




^ Reference 


35. 




^ Reference 


3e. 




° Reference 


37. 




' Cubic lattices, experimental values 


(Ref.i^). 




^ Theoretical lattice parameter, found by 


the volume optimization in the cubic 


phase, with atoms at the ideal positions. 



the cubic structure is used, a ferroelectric state is consid- 
ered using the same fractional distortions as in BaTiOa. 
Although these states may not be observed in normal 
conditions, this allows us to study the possible correla- 
tion of polarization and EFG in different systems or as a 
function of strain. 



We have considered a series of simple ABO-j, type per- 
ovskite compounds-^ For BaTiOa, PbTiOa and KNbOa 
we have considered the tetragonal experimental struc- 
tures as references. We also considered other perovskite- 
related compounds, such as BaZrOs, CaTiOs, PbZrOs, 
SrTiOs, NaNbOa and LiNbOa by considering a pseudo- 
cubic phase at the experimental lattice constants of the 
cubic paraelectric phase. The ferroelectric distortion was 
mimicked by a polar displacement of the atoms. The 
experimental displacements in tetragonal BaTiOs are 
ZT^ = 0.0203, zo\ = -0.0258, and zo2 = -0.0123, in 
fractional coordinates)^ We calculate the EFG as a func- 
tion of A, which represents the fraction of the displace- 
ments [zTii zoi, Z02) mentioned above. Therefore A = 1 
corresponds to the equilibrium (experimental) structure. 
Values of < A < 1.2 are used, i.e., A = corresponds to 
undistorted, A = 1 corresponds to the experimental equi- 
librium distortion (at 280 K) and A = 1.2 corresponds 
to 20% additional distortion. For the compounds where 



For the density functional theory calculations we used 
the projector-augmented-wave (PAW) method,— as im- 
plemented in the Vienna ab-initio simulation package 
(VASP))^ with the generalized gradient approxima- 
tion (GGA)-Perdew, Burke and Ernzerhof (PBE) func- 
tional.^^ We used a P-centered Monkhorst-Pack 7x7x7 
fc-points grid, and an energy cutoff of 400 eV. The polar- 
ization is calculated with the Berry phase approachi^ 
The EFG is calculated at the atomic sites A, B, apical 
01, and equatorial 02 sites. 

For the case of BaTiOa we also performed calcula- 
tions with the L/APW-(-lo method, implemented in the 
WIEN2K codei^ The basis set was limited by RK^ax = 
8, where RKmax stands for the product of the small- 
est atomic sphere radius Rmt times the largest iiT-vector 
Kmax of the plane wave expansion of the wave function, 
and a 6 X 6 X 6 fc-points grid was used. 



III. RESULTS AND DISCUSSION 
A. Relation between EFGs and P 

The lattice constants are shown in table HI along with 
the experimental atomic distortions in the ferroelectric 
compounds BaTiOg, KNbOg, and PbTiOg. For the other 
cases, we have considered the fractional distortion of 
BaTiOg. 

Let us start by considering the case of BaTiOs. The 
variation of the total energy with the atomic distortions 
is shown in Fig. 1(a) with the displacements considered 



along the polar axis z in both directions. The curve shows 
the expected double well profile. The energy minimum is 
reached for 85% of the experimental distortion: the slight 
discrepancy is probably due to the approximation used 
here for the exchange-correlation functional. Given that 
both the experimental/theoretical displacements are very 
small, this discrepancy is reasonable."^* In all the follow- 
ing cases, the vertical dashed and dotted lines correspond 
to the theoretical energy minimum and the experimental 
distortions, respectively. The stable state has an energy 
lower than the undistorted one by 20 meV, which com- 
pares well with previous calculations. 



Figure 1 (b) presents P as a function of the ferroelectric 



distortion. It can be seen that P is approximately a linear 
function of the distortion. This is not unexpected since 
we are considering displacive type ferroelectrics, as al- 
ready discussed in previous works.— The calculated P at 
the experimental distortion s is 2 8.6 /iCcm~^, according 
to what can be seen in Fig. 1(b) while the experimental 
one^ is 27 /zC cm~^. 

In Fig. 1(c) we show the EFG component Vzz as a 



function of distortion for all sites4i It is evident that Vzz 
shows a quadratic dependence upon the distortion. The 
sigEk^ and magnitude of the values obtained are consis- 
tent with previous calculations 4^i^ Measurements^ ob- 
tained at the O sites by NMR are in agreement for the 
equatorial site with 2.56, but at the apical site the experi- 
mental value is 2.06 while our calculated value is smaller, 
1.2 (all in units of lO^i Vm'^). 

Furthermore, we investigated how the EFG results de- 
pend on the choice of method or functional [see fig. 1(c)) . 
We have performed calculations using the L/APW-l-lo 
implementation of DFT, with the WIEN2K code. Three 
different functionals were used, which are represented 
in the plot by dotted, dashed, and dot-dashed lines, 
corresponding to the GGA-PBE^ GGA-Wu and Co- 
hen (WC),**^ and local-density approximation (LDA)^^ 
exchange-correlation functionals. The variation of EFG 
with distortion is the same, except for a small difference 
at Ti between the two implementations. There is also a 
shift in values at the O atoms, with a maximum difference 
of 0.5 X 10^^ Vm~2, when comparing LDA L/APW-l-lo 
and PBE PAW calculations, while maintaining the same 
variation with distortion. These differences are reason- 
able and the main feature, quadratic variation of Vzz with 
distortion, remains the same. All the following results are 




■2QI ^ ?? W 

-H.2 -0.8 -0.4 0.4 0.8 1.2 





-1.2 -0.8 -0.4 0.4 0.8 1.2 




500 1000 1500 

P"(10"'"C"cm"'') 

(d) 

FIG. 1. (Color online) BaTiOs in the tetragonal phase, (a) 
Energy (meV) as a function of distortion, varied in equal in- 
tervals from to 1.2. The direction of the displacements is 
shown at each side in the plot, (b) P (/xC cm^'^) as a function 
of distortion. The vertical dashed and dotted lines correspond 
to the calculated energy minimum and experimental distor- 
tion, respectively, (c) Vzz (10^^ Vm~^) as a function of dis- 
tortion (dotted, dashed and dot-dashed lines are L/APW+lo 
calculations, see text), (d) Vzz as a function of P^ (10^^^ 
C^ cm~*). Full lines are fits to the data. 
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FIG. 2. (Color online) Vzz at each site as a function of P^ 
^jQ-12 fj2 gjjj-4^ -pj^g jjjjgg j^j.g g^g ^p ^jjg Vz^P'^) data, (a) 

KNbOa, (b) PbTiOs, in the tetragonal phases. 



taken from the PAW calculations. 

In order to get rid of the distortion parameter, in 
Fig. 1(d) we plot the Vzz EFG component as a function 
of P . In this case, a fit of the data clearly shows a linear 
dependence. 

For KNbOa the trends are similar, and a linear rela- 



tion is also obtained in Fig. 2(a) In this case, the calcu- 



lated value of P at the experimental distortion is 37.4 /i 
Ccm~^, consistent with the measured value. "^^ The cal- 
culated Vzz at Nb, 2.6 x 10^^ Vm~^ is also consistent 
with 114^1 = 2.7 X 10^1 Vm"2 obtained by an NMR ex- 
periment at 220 C (Ref. Ww [considering Qsa^^ — —0.375 
(Ref.'il)]. 

For the case of FbTiOa, presented in Fig. 2(b)[ much 
larger values of polarization are obtained. The polariza- 
tion is 85 (91) at the calculated (experimental) distortion, 
higher than the measured^ 75 /i C cm^^ at room temper- 
ature. Notably large variations of the EFG in the range 
of distortions are seen at the Pb and Ti sites. The ob- 
tained values of Vzz are in reasonable agreement with Ti 
NMR experiments (Table 3 of Ref. 49) (better agreement 
is seen considering the calculation of minimum energy). 
Vzz at the 02 site has discontinuous changes, but the 
EFG tensor is continuous, as will be shown later. 

The coefficients in the Vzz{P) expression obtained are 
shown in Table HIl for all the atoms in the unit cell 
of all compounds considered. At the A, B, and api- 
cal 01 atoms, the relation found in almost all the com- 

AK,, 

pounds considered here is ( o AVy„ o ) — AP"' x 

AK^ 



( 5" 
^ 



a 



°y 



However, since these sites have rj 



0, 



the tensor is defined by only one independent parameter. 



^xx — I'l/i 



-l/2Vzz, and the quadratic coefficients also 



follow the same symmetry Qxx = o,yy = ^^/'^O'zz- The 
equatorial oxygen sites (02) do not have an n-fold rota- 
tion axis with n > 3. This implies that the asymmetry 
parameter is not zero. In this case the coefficients of P^ 
describing the variation of the EFG tensor do not show 
such a simple relation. 

The variation of the EFG with displacements can be 
understood by considering a Taylor series expansion. 



Vz, 



K°+E 



dVz: 






■Srf 



where Sri are small deviations of the atomic positions 
relative to the paraelectric structure. When the transi- 
tion involves small displacements and for atoms where 
the EFG does not undergo large changes in the transi- 
tion this expansion should converge rapidly. For sites 
without inversion symmetry, the linear term should be 
dominant, whereas for sites with inversion symmetry the 
linear term vanishes and the quadratic term in the ex- 
pansion becomes relevant;^ 

The Yzz component follows a quadratic variation for 
all the atoms in all the compounds studied, with only 
two exceptions, PbTiOa and CaTiOs. At the 02 sites 
of PbTiOa and CaTiOs there are interchanges of tensor 
components which make the description of EFG varia- 
tions in terms of the Vzz component inadequate. At the 
Pb site of PbTiOa a small quartic term in the polarization 
is found necessary for a good fit oiVzziP) {axP^+bxP'^) 
(Terms with odd powers of P are not allowed due to the 
inversion symmetry, in the paraelectric structure, of the 
sites involved in the distortion.) 

In the paraelectric structure of BaTiOa the principal 
component of the tensor Vzz for Ba and Ti is directed 
along the z axis. For the O atoms, the EFG tensor is 
also aligned with the tetragonal crystalline axes, and Vzz 
is directed to the neighboring Ti atoms. For BaTiOa, 
with increasing ferroelectric distortion (P), in all distor- 
tions calculated, the direction of Vzz remains the same 
for all atoms. The Vxx and Vyy components also main- 
tain their directions in this path, for the Ba, Ti and 01 
atoms, along the x and y axes, respectively. For the 
02 atoms, however, Vxx and Vyy do not always corre- 
spond to the same orientations. The three components 
of the EFG tensor for the 02 site in BaTiOa are shown 
At a given distortion, due to the conven- 
\Vxx\ the regular curves followed by these 



in Fig. |3(a) 
tion 



IK, 



> 



components are interchanged. For distortion A < 0.9 the 
directions for Vxx, Vyy are z and x for 02 at (0.5,0,0.5); z 
and y for 02 at (0,0.5,0.5). For higher A these directions 
are interchanged. Nevertheless, the quadratic behavior 
of EFG(P) is maintained here and this interchange is ig- 
nored. Table IIIII shows the coefficients of the quadratic 
terms for the variation of the EFG tensor components at 
each crystalline axis direction in BaTiOs. Unlike in the 



TABLE II. Coefficients a for the fits of the expression l4z(P) = V°^ + ax P^. V°^ in units of 10^^ Vm" 
VC~^ m^, for the compounds with tetragonal and cubic perovskite type structures. 



a in units of 10 





BaTiOs 


KNbOs 


PbTiOs 


BaZrOa 


CaTiOs 


PbZrOa 


SrTiOs 


NaNbOs 


LiNbOs 




^ zz 


a 


*' zz 


a 


*' ZZ 


a 


^ ZZ 


a 


K°. 


a 


v,% 


a 


^ ZZ 


a 


v?. 


a 


v°. 


a 


A 


0.16 


-1.079 


0.14 


-0.128 


5.28 


-\n%T 


0.00 


-0.945 


0.00 


-0.186 


-0.01 


-1.911 


0.00 


-0.551 


0.00 


-0.077 


0.00 


-0.012 


B 


0.78 


-0.949 


3.60 


-1.046 


5.58 


-0.562 


0.00 


-0.757 


0.01 


-0.724 


-0.02 


-0.985 


0.01 


-0.834 


0.00 


-0.163 


0.00 


-0.311 


01 


2.70 


-1.410 


2.25 


-1.042 


1.63 


-0.303 


-1.43 


-1.124 


-0.11 


-1.025 


-3.61 


-1.009 


1.08 


-1.193 


0.30 


-0.882 


0.67 


-0.940 


02 


2.32 


0.288 


1.37 


0.625 


1.18 


-0.461'' 


-1.41 


0.300 


-0.10 


0.349'= 


-3.59 


0.232 


1.11 


0.253 


0.32 


0.389 


0.69 


0.363 



'^ An additional quartic term b x P , with b — 3.626 x 10 V C^ m , is needed to get a satisfactory fit for this 

atom. 
^ This coefficient fits a component that is Vzz in the paraelectric phase, but is interchanged with other components 

in the distortion path. An additional quartic term b x P*, with b = 1.157 x 10^^ VC^^m^, is needed to get a 

satisfactory fit for this atom. 
'^ This coefficient fits a component that is Vzz in the paraelectric phase, but is interchanged with with other 



components in the distortion path. An additional quartic term b x P*, with b 
needed to get a satisfactory fit for this atom. 



-4.643 X 10^' VC-^m" 



1 


-1.31 0.019 


2 


-1.03 -0.297 


3 


2.32 0.288 



TABLE III. Fit coefficients for the 02 atoms. V{P) = 
V" + ax P^. "1/° in units of 10^^ Vm^^ a in units of 10^^ 

Compound EFG component V" a 
BaTiOs 



case of 77 = 0, the quadratic coefficients a do not have 
the same symmetry as the EFG tensor components. 

The tensor component interchanges may only happen 
at the 02 sites, where 77 changes. For BaZrOa, PbZrOs, 
NaNbOa, LiNbOa and SrTiOa there are no interchanges 
between tensor components, for KNbOs there is an in- 
terchange between Y^^ and Vyy, while for CaTiOa and 
PbTiOs there are exchanges between the three compo- 
nents. For the sake of brevity, we only show here the 02 
EFG tensor components for FbTiOa [Fig. [3(b)], where 
there are as many as five interchanges of one tensor com- 
ponent in the range of distortions considered. In par- 
ticular, there is an interchange near the experimental 
distortion of which one should be aware in EFG mea- 
surements, since a kink measured in a Vzz dependence 
could be a consequence of the convention that Vzz is the 
largest EFG component in magnitude, instead of a phase 
transition. 

Moreover, as already mentioned, to obtain a satisfac- 
tory fit to Vzz{P) at the Pb site a small quartic term 
should also be considered [Fig. 2(b)| . Vzz{P) follows ap- 
proximately a linear behavior for larger values of P^. At 
the 02 site two of the curves of components with con- 
stant direction also need quartic terms for a good fit 
[Fig. |3(b)| . The much larger displacements in this com- 
pound would indeed indicate that the Vzz{Sr) expansion 
does not converge as fast as in other cases. In order to 
confirm that with larger displacements additional terms 
must be included in Vzz{P), we performed the calcula- 
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FIG. 3. (Color online) (a) BaTiOs in the tetragonal phase: 
Vzz, Vyy, and Vxx components of the EFG tensor for the 02 
atoms as a function of P^. The lines are fits to the tensor com- 
ponents with regular variation, (b) PbTiOa in the tetragonal 
phase: Vzz, Vyy, and Vxx components of the EFG for the 
02 atoms as a function of P^. The [001] direction curves are 
also indicated. Depending on the specific atom, the other two 
curves correspond to either the [100] or [010] directions. 



tion of BaTiOa again, but this time we doubled the size 
of the distortions, allowing < A < 2.4. The variation 
of the EFG ceases to be properly described by a single 
quadratic term at the O sites, and an additional term is 
needed, as expected. 

For the other cases there is only an exception in the 02 



TABLE IV. Coefficients a for the fits of the expression V{P) — 
V° + a X P"^ for BaTiOs in the orthorhombic phase, where 
Vi,2,3 is a component of the tensor with regular variation. 
ViVa in units of 10^^ Vm-^ 01,2,3 in units of 10^^ VC^^m^ 





Ba 


Ti 


01 


02 


Vi° 


-0.17 


-0.80 


2.17 


-1.13 


ai 


0.563 


0.498 


0.318 


0.193 


Va" 


-0.11 


0.48 


-1.07 


-1.45 


02 


-0.256 


-0.829 


-0.174 


0.411 


v§ 


0.27 


0.32 


-1.09 


2.58 


as 


-0.307 


0.331 


-0.144 


-0.605 



sites of CaTiOa where a quartic term is also needed. We 
point out that CaTiOs has the smaller unit cell (TableU]), 
so the ions will be closer with the same A in comparison 
to the other compounds, which may be related to this 
exception. 

In the pseudocubic cases the chosen distortions are ar- 
bitrary. Nonetheless, we can conclude from these cases 
that, for small distortions, a simple quadratic variation of 
Vzz with polarization is seen in several different systems, 
apart from small deviations. 



B. Orthorhombic structure 
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BaTiOs exhibits monoclinic, rombohedral, and or- 
thorhombic phases at different temperatures. In order 
to see what are the differences in the EFG and P with a 
change of structure, we have made a series of calculations 
in the orthorhombic phase. We also took one experimen- 
tal measurement^ of the orthorhombic structure as the 
reference distortion (A = 1) and calculated from A = to 
A = 1.2 of this distortion keeping the lattice parameters 
(a, 6, c) constant. The theoretical distortion of mini- 
mum energy is once again found to be close to 85% of 
the experimental distortion. 

The resolved components of the EFG tensors at the 
four inequivalent atoms against P^ are displayed in 
Fig. m and the coefficients resulting from the fits are 
listed in Table IIVI In this structure there is not a single 
site with axial symmetry, and the relationship between 
EFG and P will never be as simple as discussed above. 
There are also discontinuities in Vzz at Ba and Ti sites 
due to the interchange of components. Nevertheless, ig- 
noring these interchanges, the relation for each direction 
is purely quadratic in all cases. These results, like the 
previous ones, show that there is no physical meaning in 
the change of the principal axis definition (as xx, yy, or 
zz), and that the conventional assignment may obscure 
a simpler relation with crystal axes. 
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FIG. 4. (Color online) EFG tensor components in the princi- 
pal coordinate system, of (a) Ba, (b) Ti, (c) Ol and (d) 02, 
as a function of P^ for orthorhombic BaTiOs. The lines are 
fits to the tensor components of regular variation. 



C. Correlations between EFG tensor components 

For the cases studied with tetragonal and cubic lattice 
parameters, rj j^ only for the 02 atoms. For 77 = the 
components of the EFG are trivially related, but when 
77 7^ 0, the possible correlation between tensor param- 



eters should be studied. The usual parametrization of 
Vzz and 7] assumes uncorrelated EFG components. In- 
stead, correlations between the components of the ten- 
sors may be studied in a plot of one component against 
the other. However, just plotting Vzz against Vxx results 
in a distorted cobweb plot, and EFG tensor trajecto- 
ries obtained by continuous variation of some parameter 
(temperature, distortion, pressure) may not be continu- 
ous (when 1] — 1, Vzz changes sign). Czjzek proposed a 
different system to eliminate these problems, ^"^ initially to 
deal with the analysis of amorphous systems. However, it 
is also suited to investigate the correlations between the 
tensor components t^^'^^ This plot uses —2Vxx as a func- 
tion of \2{2Vzz + Vxx)/V3\- With this linear combination 
of tensor components, the tractories are always continu- 
ous, and are straight lines if there is a linear dependence 
oiVzz on Vxx- 

The trajectories in the Czjzek plot of all the tetrag- 
onal and pseudocubic cases are shown in Fig. [5] for the 
equatorial oxygen atoms, where i] changes and a non- 
trivial correlation may be present. In this plot, the lines 
of constant 77 are the lines emerging from the origin: the 
boundary lines correspond to 77 = and the horizontal 
line corresponds to 77 = 1. The herringbone lines corre- 
spond to constant Vzz ■ Reflections at the boundary of the 
plots, seen for example in BaTiOa and KNbOa, are asso- 
ciated with the interchange of Vxx and Vyy components 
previously shown, and the crossing of the trajectory of 
the line 77 = 1 corresponds to a change in sign and orien- 
tation of Vzz ■ The equatorial oxygen sites still have axial 
symmetry (rj — 0) in the paraelectric structure in the cu- 
bic structures, since there is four-fold rotation symmetry 
around the axes connecting the O to the B sites. There- 
fore, in these cases, the trajectories start at the boundary 
of the plot. For the tetragonal cases, this symmetry is 
lost and i] ^ even in the paraelectric phase: rj — 0.12 for 
BaTiOg, Tj = 0.33 for KNbOg and 77 = 0.93 for PbTiOa. 
The upper (lower) wedge is for positive (negative) Vzz- 
The values of Vzz{02) change significantly for the dif- 
ferent compounds, even changing sign. KNbOa, SrTiOa, 
LiNbOa and NaNbOa always have positive values of Vzz , 
while the zirconates BaZrOs and PbZrOa have negative 
values of Vzz in distortions considered. The EFG changes 
sign in the trajectories of CaTiOa and PbTiOa. In the 
case of CaTiOa the trajectory overlaps with itself. It 
starts at negative EFG and 77 = 0, then it goes to pos- 
itive EFGs, in the direction of the arrow, it is reflected 
and comes back the same way, changing sign again, and 
finally it is reflected to the horizontal path. In general 
the positive EFGs increase in absolute value with dis- 
tortion, while the negative ones decrease. Exceptions in 
this respect are parts of the trajectories in PbTiOa and 
CaTiOa. 

The lengths of the trajectories in this plot are pre- 
sented in Table |Vl The lengths are markedly larger for 
PbTiOa, KNbOa, LiNbOa, and NaNbOg. The cases of 
PbTiOa and KNbOa correspond to different fractional 
distortions in comparison to the other cases. However, 



TABLE V. Approximate values for the lengths of the 
EFG(02) trajectories in the Cjzjek plots, relative to the 
length of BaTiOa , and their asymptotic asymmetry parame- 
ter (r/oo). 



Compound length rjoo 

BaTiOs 1 0.92 

SrTiOs 0.99 1 

CaTiOs 0.85 1 

BaZrOs 0.88 0.81 

LiNbOa 1.91 0.63 

NaNbOs 2.02 0.77 

PbZrOs 0.59 0.92 

KNbOa 2.87 1 

PbTiOa 5.76 0.46 



all the other cases have fractional distortions equal to 
BaTiOa. In this respect, it is interesting to remark that 
both niobates (LiNbOa and NaNbOa) have trajectories 
with approximately twice the length of the others, in- 
dicating that the local charge at the 02 atoms is more 
sensitive to displacements when Nb is the B site. 

From these plots it is also of interest to determine the 
"asymptotic" asymmetry parameter (rjoo)^ related to 
the slope of the trajectories after reflections at the bound- 
aries, corresponding to the limiting value of 77 for a given 
trajectory. These values are presented in Table |Vl It can 
be seen that 7700 is large in most cases, with a smaller 
value for PbTiOa. It should be related to the different 
features of the distortion, due to the fact that this com- 
pound has ferroelectricity driven by the lone pair of s 
electrons at the Pb ions, while ferroelectricity in other 
compounds is related to the (P configuration at the B 
site. 

If the trajectories in this plot are straight lines, this 
shows that the EFG tensor components are linearly re- 
lated. This is usually the case, except for LiNbOa and 
NaNbOa, where the trajectory is approximately straight 
for low values of P, but becomes curved for higher values. 
The linearity in almost all cases means that the whole 
tensor, when considering the distortions that give rise to 
ferroelectric polarization, can be described by a single 
parameter. This might have important implications on 
the way experimental data should be analyzed, since the 
components are not independent, and the combination 
of all information will be more useful than each of them 
taken separately. The study of the EFG dependence on 
temperature (or on other variables) usually done by sep- 
arating Vzz and ?7, can be performed by identifying the 
global single parameter. 
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FIG. 5. (Color online) Cjzjek plots of the EFG tensor com- 
ponents at the equatorial 02 sites with P (or distortion) as 
the implicit parameter in the perovskite materials studied, 
covering < A < 1.2, in units of 10^^ Vm~^. The arrows 
indicate the "initial" point and direction of the trajectories, 
corresponding to P = 0. Trajectories start at r; = except 
for the tetragonal cases. 
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FIG. 6. (Color online) (a) T4z at the A site as a function of 
P^, showing the different coefficients (slopes of straight lines), 
(b) Viz as a function of P^ for PbTiOs. (c) Coefficients of 
quadratic term of Vzz(P), at the A site, for a series of ma- 
terials with the perovskite structure, as a function of atomic 
number Za- The line is a quadratic fit. 



IV. RELATION BETWEEN EFG AND P IN 
DIFFERENT MATERIALS 



In the following we analyze the relations found in differ- 
ent materials. The quadratic coefficients a in the Vzz{P) 
expression follow an interesting trend with atomic num- 
ber at the A sites (Za), shown in Fig. |6l Figures 6(a) 
and |6(b)| show Vzz at the A site for the different cubic and 
tetragonal compounds. PbTiOa introduces larger varia- 
tions for P^ and Vzz which may follow from the larger 
displacements in the experimental distortion, related, as 
previously discussed, to the different nature of the ferro- 
electricity. 

Vzz is for the paraelectric cubic cases, while for the 



paraelectric tetragonal structures Vf^ has small nonzero 
values, except PbTiOs. The slopes of the trajectories cor- 
respond to the coefficients of the quadratic term, which 
are shown in Fig. 



6(c; 



The A sites with a larger num- 
ber of electrons show a larger magnitude of a, that fol- 
lows a quadratic dependence on Za- There are two pairs 
of compounds with the same Za, but different metal 
B sites, namely BaTiOs/BaZrOs, and PbTiOa/PbZrOs. 
The same trend is still qualitatively followed although 
small changes in the values of the coefficient are seen, 
and the changes do not follow a pattern: a is higher for 
B = Ti than B ~ Pb in the cases with A — Ba, but the 
opposite happens when A = Pb. We did not find other 
trends for the other coefficients, or for V^^, as a function 
of the atomic number of the different sites. 



CONCLUSIONS 
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The EFG, working as a local analog of P, has the 
added advantage of higher spatial resolution and short 
time scales, which allows local probing of nanoscale phe- 
nomena at specific lattice sites. The different types of lat- 
tice sites, including defects, can be discriminated by their 
different EFGs. Its critical behavior in phase transitions 
may be analyzed, and the order of the transitions can 
be established with a high degree of detail. In hysteresis 
loops, it can act as a static measurement of the electric 
polarization in individual domains, not limited by depo- 
larization effects and is much more sensitive to polariza- 
tion reversal. It is also suited to probe phase coexistence 
or inhomogeneous polarization states with atomic selec- 
tivity, well beyond the reach of conventional polarization 
measurements. Moreover, piezoelectric force microscopy 
is usually restricted to studies near the surface of sam- 
ples, while EFG studies can be performed in bulk or at 
the surface of materials, by using diffusion, evaporation 
or implantation techniques to add the probe atoms in the 
environments to study. 



Therefore, we hope this work will stimulate more EFG 
studies. However, one limitation of our results should be 
mentioned. In this paper we have analyzed the variation 
of the EFG with the scaling of the ferroelectric distor- 
tions as a whole. However, in experiments, the parame- 
ters, such as temperature, in general produce variations 
involving, for example, additional distortion modes, lat- 
tice vibrations, and lattice expansions. The variation of 
each independent atomic displacement may have a com- 
plex variation with changing temperature. Therefore, the 
relation investigated here can not be directly used to infer 
polarization variations from EFGs or vice-versa, except 
in well understood cases. Further work to improve this 
limitation could involve the study of different structural 
changes or temperature effects. 
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